Kramer-Pesch approximation for analyzing field-angle-resolved measurements made in 
unconventional superconductors: A calculation of the zero-energy density of states 
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By measuring angular-oscillation behavior of the heat capacity with respect to the applied field 
direction, one can detect the details of the gap structure. We introduce the Kramer-Pesch approx- 
imation (KPA) as a new method to analyze the field-angle-dependent experiments quantitatively. 
We calculate the zero energy density of states for various combinations of typical Fermi surfaces 
and superconducting gaps. The KPA yields a merit that one can quantitatively compare theoretical 
calculations with experimental results without involving heavy numerical computations, even for 
complicated Fermi surfaces. We show an inadequacy of the frequently-used Doppler-shift technique, 
which is remedied by application of the KPA. 

PACS numbers: 74.20.Rp, 74.25.Op, 74.25.Bt 



Much attention has been paid to investigation on 
the mechanism of exotic superconductivity emerged in 
highly-correlated materials like high-T c cuprates and 
heavy-fermion compounds, in metallic materials like non- 
magnetic borocarbides, and in organic materials like 
BEDT-TTF salts. Many exotic superconductors (SCs) 
have zeros of the single-particle excitation gap called 
nodes on the Fermi surface (FS). It is then important 
to identify the detailed gap structure; either line or point 
nodes and their locations on the FS. The identification 
of them is crucial for understanding the unconventional 
pairing mechanism of exotic SCs [J 

To this end, various experimental methods have been 
developed. The power-law temperature dependence of 
physical quantities such as specific heat and nuclear mag- 
netic relaxation rate provides partial information on the 
gap structure, namely a difference between line and point 
nodes. In the last decade, the angular dependence on an 
applied magnetic field has facilitated us deducing more 
directly the gap structure including the information on 
the positions of gap nodes 0, One of such power- 
ful methods is the angular-resolved heat capacity mea- 
surement . By measuring angular-oscillation behavior 
of the heat capacity C(H) with respect to the applied 
field H direction, one can detect the details of the gap 
structure, especially for the locations of its nodes. The 
angular-resolved thermal conductivity is also such an ex- 
perimental technique These measurements are capa- 
ble of achiving high precision enough to discuss details of 
the oscillatory structure. Owing to potential abilities to 
reveal the gap structure in much more various exotic SCs, 
experimental and theoretical progresses in this direction 
are highly influential developments in the condensed mat- 
ter physics. 

For instance, it is believed that YM2B2C is the first 
important material in which superconductivity with a 
point- node gap was confirmed experimentally 0, [H] . One 
of the evidences of point nodes in this material has been 



considered to be an experimental result that the polar- 
angle dependence of the thermal conductivity is strong 
as shown by Izawa et al. in Fig. 2 of Ref. [4[ (and in 
Fig. 7 of Ref. (H). They proposed that the polar-angle 
dependence would be weak if YM2B2C had line nodes, 
instead of point nodes, as shown in Fig. 4(b) of Ref. [3[. 
Their proposition is based on a Doppler-shift (DS) anal- 
ysis 0, [f for an ellipsoidal FS @. The DS method, 
which was crucial for that proposition of point nodes in 
YM2B2C, has been frequently utilized to analyze data 
also for other various materials [3J] . As noted below, how- 
ever, the DS method faces certain difficulties, making 
it highly questionable for quantitative analyses, unfortu- 
nately. It should also be noted that YNi2B2C has highly- 
anisotropic complicated FSs Q in comparison to rather 
simple models of FSs adopted in usual studies. 

There is a twofold problem in the DS method for an- 
alyzing experimental results. One of the problems is an 
overcstimation of the angular-oscillation amplitude. The 
DS method gives rise to an inflated oscillation amplitude 
of the order 30 % of the maximum amplitude [8[, al- 
beit actual oscillation amplitudes are only of the order a 
few percents in many experiments (2j. Indeed, a full nu- 
merical calculation [9( leads to an oscillation amplitude 
consistent with those experiments, suggesting an inade- 
quacy of the DS method. The other problem is related to 
an issue whether or not a cusp-like minimum appears for 
the magnetic field H || gap-nodal direction. Miranovic et 
al. [9( showed that a broad (non-cusp-like) minimum ap- 
pears in a full numerical calculation for a spherical FS. In 
experiments, the cusp-like minima do not appear in most 
materials j2j. Only in YM2B2C the cusps were observed 
[H As discussed later, for addressing this issue the 
DS method appears to be inadequate in materials with 
complicated FSs because of an approximation inherent 
in its formulation. Although full numerical calculations 
could resolve both of these problems in principle, they re- 
quire significantly long computational time and it is hard 
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to conduct such direct calculations in systems with re- 
alistically complicated FSs. The Brandt-Pesch-Tewordt 
(BPT) method is alternatively proposed, which is an ap- 
proximation adequate near the upper critical fields H C 2 
[10(. The BPT is also not an easy-to-calculate method 
in low magnetic fields H and at low temperatures T, 
unfortunately. In this respect the DS method has been 
a favorably easy tool, but it faces the above mentioned 
problems. Therefore, it is now highly desirable to develop 
a new tractable method that provides quantitatively cor- 
rect results. 

In this Letter, we first introduce the Kramer-Pesch ap- 
proximation (KPA) as an efficiently applicable method to 
analyze the angular- resolved experiments. We will show 
that the KPA yields the zero energy density of states 
(ZEDOS) consistent quantitatively with results of direct 
numerical calculations. The computational time required 
for the KPA is almost the same as for the DS method 
and is significantly less than for direct numerical calcula- 
tions. Using the KPA, one can also obtain the density of 
states (DOS) in a system with realistically complicated 
FSs without involving heavy numerical computations. In 
addition, we will show that the DS method is inadequate 
and the KPA would be indispensable for analyzing ex- 
perimental results, since some information on the FS is 
neglected in the DS method. 

The DOS is the basis for analyzing physical quanti- 
ties such as the specific heat and the thermal conductiv- 
ity (Tfj| . It is obtained from the regular Green function 
g within the quasiclassical theory of superconductivity, 
which is represented by a parametrization with a and 
b as g = —(1 — a6)/(l + ab). They follow the Riccati 
equations (h = 1) [111 ]: 

v F ■ Va + 2Q n a + aA*a - A = 0, (1) 
v F ■ Vb - 2u n b - bAb + A* = 0. (2) 

Here, vf is the Fermi velocity, and iu> n = iw n + (e/c)vF-A 
with the Matsubara frequency uj n and the vector poten- 
tial A. 

First, we show what is the difference between the 
DS method and the KPA from a fundamental view- 
point. In the DS method, we separate the pair po- 
tential A into the amplitude and the phase &(x,k) as 
A(x,k) — |A|exp(i$), and consider that the solution 
of Eq. ([T]) is written as a — exp(— i&)a 12]. Here, x 
is the coordinate along the direction of v-p. Assuming 
dd/dx — 0, the equation can be exactly solved in an 
analytical way and the solution is equivalent to that in 
the bulk region replacing the energy uj n by the Doppler 
shifted energy oj n + ivp ■ mv s . The DS method is based 
on the fact that the energy of a quasiparticle (QP) shifts 
by the Doppler-shift energy 8E = mv-p ■ v s due to circu- 
lating supercurrents [l3j , where v s means the superfluid 
velocity around a vortex and v s _L H . Consequently, 
the ZEDOS is obtained as ~ Re(\SE\/ ^{SE) 2 - |A| 2 ) 
[1, [13, [3], where the angular brackets mean averaging 
over both a FS and a unit cell of vortex lattice. For H | 
gap-nodal direction, near gap nodes the QPs whose Fermi 



velocities are parallel to H (vf || H) lead to SE = 0. 
As a result, those QPs are excluded from the contribu- 
tion to the ZEDOS. Therefore, C(H)/T, which is pro- 
portional to the ZEDOS in the limit T — > 0, has a min- 
imum for H || gap- nodal direction |8|. The DS method 
(da/dx = 0) is an approximation neglecting the spatial 
variation of \a\ = \a\, so that it breaks down near a vortex 
core, around which the pair potential A, a, and b sub- 
stantially change. The KPA avoids this neglect. Indeed, 
Dahm et al. [12[ showed that the DS method misses im- 
portant contributions from vortex-core states extending 
into the gap-nodal directions far away from a vortex core. 

To take account of those contributions, in the KPA we 
expand the Riccati equations up to first order in the im- 
pact parameter y and the energy u> n (ioj n — > E + irj and y 
is the coordinate along the direction perpendicular to vf) 
[IH [I(| . It follows that only the Andreev bound states 
around a vortex core is taken into account by the KPA 
and this approximation is appropriate in the low-energy 
region \u> n \ <C |A|. The KPA is based on the fact that the 
delocalized states due to the QPs near gap nodes can be 
regarded as the Andreev bound states spreading to infi- 
nite distances, so that the Andreev bound states predom- 
inantly contribute to the DOS in the low energy. Indeed, 
Kopnin and Volovik [k| confirmed that, for the DOS, 
the DS method and a calculation by considering only the 
bound states (corresponding to the KPA here) yielded 
results of the same order in the case of two-dimensional 
(2D) isotropic FS under a magnetic field applied perpen- 
dicular to the 2D plane. It should also be noted that 
the KPA works well widely far away from a vortex core 
since the relation \uj n \ <C \A\ is marginally satisfied even 
near gap nodes in the low energy and the KPA correctly 
describes vortex core states extending into the gap-nodal 
directions [l6[. Hence, the KPA appears to be adequate 
for calculating the low-energy DOS in the vortex state. 

By the expansion procedure based on the KPA [l6j| , an 
analytic solution of the Riccati equations can be obtained 
around a vortex. The DOS is calculated as (— Reg R ) with 
the retarded Green function g R accordingly. We find the 
new expression for the angular-resolved DOS in general 
form, 




dS F A[cosh(ir/£o)] 



2tH£ \J \vf\ {E-E y y-t V - / g] 

.(3) 

Here, the azimuthal (polar) angle of H is «m (0m) m a 
spherical coordinate frame fixed to crystal axes, and cISf 
is an area element on the FS [e.g., gLSf = k F sin 8d(f>d8 
for a spherical FS in the spherical coordinates (k,cf),8), 
and dSF = kF a bd(f>dk c for a cylindrical FS in the cylindri- 
cal coordinates (k a b, 4>, k c )]. In the cylindrical coordinate 
frame (r, a, z) with z || H in the real space, the pair po- 
tential is A = AoA(A;f) tanh(r/£o) exp(ia) around a vor- 
tex, Ao is the maximum pair amplitude in the bulk, A — 
A|, and (--^sp = J r ° rdr ^ ■ ■ ■ da / (trr 2 ) is the real- 
space average around a vortex, where r a /£o = y/H^/H 
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Line-nodes with the spherical FS 



(a) KPA calculation 



(b) Numerical and KPA calculations 
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FIG. 1: (Color online) Angular dependence of the ZEDOS 
calculated by the KPA in a d x 2_ y 2-w&ve SC with the spherical 
FS. Circle o, box □, and triangle A indicate the results for 
the polar angles 8m = 7r/2, 7t/3, and 7r/4, respectively. The 
smearing factor r\ = 0.05 An and the cut-off length r a = 7£o- 
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FIG. 2: (Color online) Angular dependence of the ZEDOS 
obtained by (a) the KPA and (b) the direct numerical and 
KPA calculations in the d a .2_ H 2-wave SC with the cylindrical 
FS. (a) The magnetic field tilts from the c axis by polar angle 
0m- V = 0.05A and r a = 7£ - (b) 6 M = tt/2, r) = 0.001A , 
and r a = 5£o- 



[H c2 EE $o/(^ 2 ), $0 = 7TT»fl]. 



rcos(a — 9 V ), 



y = rsm(a~9 v ), and E y = A X 2 y/(^ h). 8 v (k F ,a M ,0 M ) 
is the angle of v-p± in the plane of z = 0, where a and 
9 V are measured from a common axis [lil . fl6| . i>f_l 
is the vector component of vp(kp) projected onto the 
plane normal to H — (aM,fti). |fF±(fep, QiMi #m)| = 
fFo(aM^M)^(fcF. aM,^) and t>Fo is the FS average of 
I^fjJ [14| - £o is defined as £o = vfo/(ttA ). We consider 
here a clean SC in the type-II limit. The impurity ef- 
fect can be incorporated through the smearing factor rj. 
We should note that the expression by the KPA includes 
only the contribution of the Andreev states as mentioned 
above and it excludes the bulk DOS under the zero mag- 
netic field. This is, however, sufficient for considering 
the H dependence of the DOS in low energies E because 
the bulk DOS does not depend on H and the Andreev 
states dominate in low E. While Eq. © might seem 
more complicated than the DS expression, a load of re- 
quired numerical integrations is of the same order and is 
much less than direct numerical calculations of the Ric- 
cati equations efficiently. 

We show results for two cases comparing the KPA 
with direct numerical calculations to prove the quanti- 
tative correctness of the KPA. We consider a d x 2_ y 2- 
wave SC. First, we calculate the ZEDOS (E = 0) by 
the KPA for the spherical FS (and the ri-wave pairing 
A = cos 2cf) sin 2 0). As shown in Fig. [TJ the cusp- like min- 
imum does not appear and the oscillation amplitude is of 
the order 3.5 %, coinciding well with the full numerical 
results by Miranovic et al. 0] where no cusp-like mini- 
mum appears and the amplitude is of the order 3.5 %. 
We have also conducted a DS calculation under the same 
condition and obtained a result that there is indeed no 
cusp, but the amplitude is of the order 4.9 % with less 
coincidence. Second, we calculate the ZEDOS by the 
KPA for the cylindrical FS (and A = cos 2(f)). As shown 
in Figs. Ufa) and[2]Jb), the cusp- like minimum appears 
and the oscillation amplitude is of the order 20-23 %. 
We also calculate the ZEDOS numerically by solving di- 
rectly the Riccati equations to compare with the KPA 



(a) Point-nodes with the cylindrical FS (b) Point-nodes with the spherical FS 
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FIG. 3: (Color online) Angular dependence of the ZEDOS 
calculated by the KPA in SCs with a point-node gap on (a) 
the cylindrical and (b) the spherical FS. r\ = 0.05An and 
r a = 7£o- 



result, and find that the cusp structure appears and the 
oscillation amplitude is of the order 18 % as shown in 
Fig. HJb). The DS method, in contrast , gives rise to a 
too large amplitude of the order 30 % [8j. In summary 
the above two cases show that the results of the KPA is 
quantitatively equivalent to those of the direct numeri- 
cal calculations. We have also confirmed that within the 
KPA the dependence of the ZEDOS on the field strength 
shows H 5 behavior [13j in gap-nodal SCs. 

We consider the polar-angle dependence of the oscil- 
lation amplitude of the ZEDOS, motivated by the ex- 
periments on YM2B2C. We discuss the difference of the 
shape of FSs and of the gap structure. We calculate the 
ZEDOS for the following four combinations; SCs with 
the line-node or the point-node gap on the cylindrical or 
the spherical FS. Here, the point-node gap is the s + g 
wave 3 A = (1 + sin 4 f3cos40)/2 for the spherical FS 
and is A = (I + cos 4 (fc z /2) cos4</>)/2 for the cylindrical 
FS. As shown in Figs. Q] and [31(b) , for the spherical FS the 
polar-angle dependence in the line-node case [Fig. Q] is 
different from that in the point-node case [Fig.[3Jb)]. The 
polar-angle dependence is relatively weak in the case of 
the line- node gap and spherical FS [Fig. [T]. In contrast, 
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FIG. 4: Schematic figures of (a) the spherical and (b) the 
cylindrical FS in the d-wave case. The left and right panel 
denote a top and a side view, respectively. The hatched re- 
gions indicate the low-energy QPs around gap nodes. 

as shown in Figs.^a) and[3Ja), for the cylindrical FS the 
polar-angle dependence is strong both in the cases of line 
nodes and point nodes. These results suggest, that one 
cannot readily determine whether a SC has point nodes 
or line nodes on the basis of the polar-angle dependence 
of angular-resolved measurements, since the polar-angle 
dependence of the ZEDOS depends more eminently on 
the shape of FSs than on the type of gap nodes. 

We understand the difference in the results between 
the spherical and cylindrical FSs in the d-wave case as 
follows. The QPs near gap nodes (the nodal QPs) pre- 
dominantly contribute to the low-energy DOS. The QPs 
with Dp parallel to H are excluded from the contribution 
to the DOS for the same reason as mentioned before, so 
that the ZEDOS becomes small when Dp of the nodal 
QPs becomes parallel to H . Now, on the spherical FS, 
the nodal QPs with vp parallel to H are located only 
in narrow pole regions, and therefore the decrease of the 
ZEDOS due to such a QP exclusion is not so significant. 
In contrast, on the cylindrical FS (see, Fig. BJb)) the 
nodal QPs with v-p \\ H are located around long line-like 
regions. The QP exclusion occurs in such wide regions, 
so that the ZEDOS dramatically decreases for H || gap- 
nodal direction. As a result, the ZEDOS prominently 
changes near this gap-nodal direction and such behavior 
leads to the cusp-like structure as in Fig. [2l One can 
understand the difference in the polar-angle dependence 
in the same way. 

We show that the DS method is inadequate to analyze 
the angular-resolved experiments. In the DS method, 
only the nodal QPs that satisfy the relation \mvp -v s \ > 
A | contribute to the ZEDOS as noticed from the DS ex- 
pression. That is, the contribution of the anti-nodal QPs 
in the region where \mv-p ■ v s \ < A is completely ne- 
glected. This means that the DS method overestimates 
the relative contribution of the QPs near gap nodes on 
the FS. As a result, the contrast in the angular- resolved 
ZEDOS between the nodal and anti-nodal directions is 
exaggerated, leading to too large oscillation amplitudes 
as mentioned before. This problem in the DS method 
would be fatal when analyzing the angular-resolved ZE- 
DOS for general SCs with complicated FSs where the 
contributions of various FS regions are equally expected 



to be important. For example, if a FS has a narrow re- 
gion where its shape is similar to the cylindrical FS, a 
point-node gap in this region gives rise to unexpected 
cusp structure similar to the cusped ZEDOS for a line- 
node gap because the DS method considers only regions 
near gap nodes and does not clearly distinguish whether 
those gap nodes are distributed in side-by-side array (i.e., 
a line node) or in isolation (i.e., point nodes). Indeed, we 
have confirmed that a DS calculation under the same con- 
dition as in Fig. EKa) gives rise to inappropriate cusp-like 
structure. In contrast, the whole region on a FS is consid- 
ered by the KPA because, through the vortex-core states, 
the KPA takes account also of the QPs in the anti-nodal 
region on the FS, which are neglected in the DS method. 
In other words, the DOS obtained by the KPA includes 
correctly both the contributions of the "gap node" in the 
real space (i.e., the vortex core) and the gap nodes on a 
FS in the momentum space. 

The cusp-like minimum was observed in YM2B2C 
[3, 0- From the results for the typical FSs shown in 
Figs. it is anticipated that, as realized in the case 
of the cylindrical FS with the line-node gap, the cusp- 
like minimum appears in the case when there exist many 
nodal QPs with vp directed in a common direction. This 
suggests that, if YM2B2C has a point-node gap, the di- 
rectional distribution of Vf around gap nodes should have 
singular structure under favor of the highly anisotropic 
FSs [7| in this material. This anticipation would be con- 
sistent with a claim [l4j that the point nodes would be sit- 
uated in the regions on the FS connected by nesting vec- 
tors. The importance of the FS structure has been em- 
phasized by Udagawa et al. with the BPT calculation [l7j. 
Further detailed analyses based on the realistic FSs are 
necessary to gain conclusive insight, for which the main 
result of the present work, Eq. ©, applicable to general 
shapes of FSs will play an efficient role. 

In conclusion, we introduced the KPA as a new method 
that allows us to analyze angular-resolved experiments 
without heavy numerics. According to the results for the 
typical FSs, the FS shape eminently affects the angular- 
resolved ZEDOS, suggesting that one cannot readily 
identify a superconducting gap in a material without 
considering its realistic FS shape. The KPA is a use- 
ful method as a more proper approximation in a low-H 
and low-T region than the DS technique and is comple- 
mentary to the BPT method appropriate near H C 2. 
We expect that the present work will stimulate further 
experimental efforts to detect precise oscillatory struc- 
ture. 

We thank Y. Kato, K. Izawa, M. Udagawa, P. Mira- 
novic, N. Nakai, M. Ichioka, K. Machida, and M. Machida 
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